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ABSTRACT 

The spherical Mexican Hat wavelet is introduced in this paper, with the aim of 
testing the Gaussianity of the Cosmic Microwave Background temperature fluctua- 
tions. Using the information given by the wavelet coefhcients at several scales, we 
have performed several statistical tests on the COBE-DMR maps to search for ev- 
idence of non-Gaussianity. Skewness, kurtosis, scale-scale correlations (for two and 
three scales) and Kolmogorov-Smirnov tests indicate that the COBE-DMR data are 
consistent with a Gaussian distribution. We have extended the analysis to compare 
temperature values provided by COBE-DMR data with distributions (obtained from 
Gaussian simulations) at each pixel and at each scale. The number of pixels with 
temperature values outside the 95% and the 99% is consistent with that obtained 
from Gaussian simulations, at all scales. Moreover, the extrema values for COBE- 
DMR data (maximum and minimum temperatures in the map) are also consistent 
with those obtained from Gaussian simulations. 



1 INTRODUCTION 

Cosmic Microwave Background (CMB) temperature fluctu- 
ations are predicted to be Gaussian by standard Inflation- 
ary theories. Topological Defect models and non-standard 
Inflation can produce non-Gaussian CMB temperature fluc- 
tuations. The most recent results of CMB observations by 
Boomerang, DASI and MAXIMA-I (Netterfield et al. 2001, 
Pryke et al. 2001, Stompor et al. 2001) seem to indicate that 
any possible non-Gaussianity present in the data will more 
likely be produced by non-standard Inflationary models. De- 
termination of the possible degree of non-Gaussianity of the 
temperature fluctuations in CMB maps is therefore a funda- 
mental question to be answered. Several tests have been pro- 
posed on real space ranging from calculating skewness and 
kurtosis (Luo & Schramm 1993a), 3-point correlation func- 
tions (Luo & Schramm 1993b, Kogut et al 1996), Minkowski 
functionals (Gott et al. 1990, Torres et al. 1995, Kogut et al. 
1996, Schmalzing & Gorski 1997, Novikov et al. 2000) and 
topological quantities related to characteristics of excursion 
sets above a given threshold (Bond & Efstathiou 1987, Vit- 
torio & Juzkiewicz 1987, Martinez-Gonzalez & Sanz 1989, 
Barreiro et al. 1997,1998,1999), to multifractal analysis, par- 
tition function (Pompilio et al. 1995, Diego et al. 1999) and 



surface rougness based analysis (MoUerach et al. 1999). Us- 
ing the expansion of temperature fluctuations in spherical 
harmonics, one can test the non-Gaussianity of the CMB 
by calculating the bispectrum (Luo 1994, Heavens 1998, 
Ferreira et al. 1998, Magueijo 2000). Moreover, several re- 
cent papers have appeared proposing analysis in wavelet 
space (Pando, Valls-Gabaud & Fang 1998, Hobson, Jones 
& Lasenby 1999, Barreiro et al. 2000, Aghanim, Forni & 
Bouchet 2000). 

Many of the works cited above aimed to detect non- 
Gaussianity in the 4-year COBE-DMR data, which currently 
remains the only publicly available full-sky CMB data set. 
Ferreira, Magueijo & Gorski (1998), using a reduced bis- 
pectrum that did not consider correlations between differ- 
ent scales, detected a non-Gaussian signal which was most 
likely produced by a systematic artifact in the COBE-DMR 
data (Banday, Zaroubi & Gorski 2000). No detection of 
non-Gaussianity was obtained after that artifact was re- 
moved. However, non-Gaussianity has been detected at the 
97% — 99.8% level by using an extended bispectrum that 
takes into account correlations between different scales by 
Magueijo (2000). 

Non-Gaussianity has not been revealed by any of the 
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other tests applied by different authors, including those car- 
ried out in wavelet space (Mukherjee, Hobson & Lasenby 
2000, Aghanim, Forni & Bouchet 2000). Although the 
COBE-DMR data cover the whole celestial sphere, most 
of these works studied the wavelet coefficient distributions 
corresponding to pixels covering Faces and 5 in the Quad- 
Cube pixelisation (White & Stemwedel 1992) . Therefore 
the wavelet coefficients are calculated on the plane. The 
first attempt to detect non-Gaussianity in the COBE-DMR 
maps analysing all the available pixels and therefore using 
a wavelet projected on the surface of the sphere was pre- 
sented by Barreiro et al. (2000). The wavelet used in that 
work was the spherical Haar wavelet (as introduced to CMB 
analysis by Tenorio et al. 1999). Skewness, kurtosis and two- 
scale correlation functions were computed in wavelet space 
corresponding to different scales for the 4-year COBE-DMR 
data. The tests could not confirm any strong evidence of 
non- Gaussianity. 

Our aim in this paper is to attempt to detect non- 
Gaussianity in the 4-year COBE-DMR data by applying 
different statistical tests to the wavelet coefficients obtained 
by convolving the COBE-DMR map with a stereographic 
projection of the Mexican Hat wavelet on the sphere. This 
is the first time that the Mexican Hat wavelet has been used 
in this projection. The definition and basic properties of this 
wavelet, in particular in relation to its application on the 
sphere, are presented in Section 2 (an extended discussion 
will be given in Martmez-Gonzalez et al. 2001). 

It should be noted that, while the Haar wavelet privi- 
leges some directions (with diagonal, vertical and horizontal 
components), the Mexican Hat wavelet is isotropic. We will 
compare the ability of both wavelets to detect non-Gaussian 
features in the COBE-DMR data. The method and tests ap- 
plied in this paper are presented in Section 3. Section 4 is 
dedicated to present and discuss the results. Conclusions are 
also included in that section. 



2 THE SPHERICAL MEXICAN HAT 
WAVELET 

In many cases of astrophysical interest, data are given on 
the sphere. In this paper, we are interested in the analy- 
sis of the COBE-DMR data in the HEALPix pixehsationQ 
(Gorski, Hivon & Wandelt 1999), but current and future 
experiments measuring CMB anisotropy can also be inves- 
tigated using these methods. A conventional analysis tech- 
nique is based on spherical harmonic decomposition of the 
temperature anisotropy for which local information is lost. 
On the other hand, wavelets -defined on the line/plane- have 
been extensively used (spectral and image analysis) in astro- 
physical applications during the last years. In particular, the 
continous isotropic wavelet transform of a 2D signal f{x) is 
defined as 



w{R,b)= I dxf{x)^{R,h;x), 



being 
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and where w{R,b) is the wavelet coefficient associated to 
the scale R at the point with coordinates b. The limits in 
the double integral are -co and oo for the two variables, tp 
is the wavelet mother function that satisfies the constraints 
J dxjp — Ojjdxip'^ — 1, and the admissibility condition 
= (27r)^/ dkk~'^ tp'^{k) < oo, where Tp{k),k = \k\, is the 
Fourier transform of ip{x). A reconstruction of the image 
can be achieved with the inversion formula 



1 



dRdb^w{R,b)i^i^-^-^). 



(3) 



The analysis of 2D images can be done with the so- 
called Mexican Hat wavelet (MEXHAT) . It is given by 



■ip{x) = 



(27r) 



1/2 



(2 



V2 



(4) 



that is proportional to the laplacian of a Gaussian. It has 
been sucessfuUy applied to extract point sources from CMB 
maps (Cayon et al. 2000, Vielva et al. 2000). 

The problem arises when one tries to extend the conti- 
nous isotropic wavelet to the sphere. There had been some 
proposals to do such an extension (e. g. analyzing functions 
defined through spherical harmonics, Freeden 1997; second 
generation wavelets, Sweldens 1996). However, one would 
like to incorporate the following properties: i) the trans- 
form is defined by a wavelet, i. e. J dfl^s^ = 0, ii) the 
wavelet transform must involve translations and dilations 
on the sphere and iii) for small scales the transform must 
reduce to the usual continous isotropic wavelet on the plane. 
Recently, Antoine & Vandergheynst (1998) have presented a 
group-theoretical derivation of the continous wavelet trans- 
form on the 2-sphere S2 satisfyng the three previous con- 
ditions. They identify the appopriate transformations as: i) 
translations are given by elements of the 3D rotation group 
5*0(3); ii) dilations around any point are the inverse stereo- 
graphic projection of standard dilations on the tangent plane 
to the point. Using polar coordinates (6^, (j)) around any fixed 
-but arbitrary- point on 5*2, such a projection is defined by 



x = 2 tan - 



(5) 



where {x, cj>) are polar coordinates on the tangent plane. 
Moreover, the analyzing wavelet, with respect to any fixed 
but arbitrary point, is given by (Martmez-Gonzalez et al. 
2001) 



*s,(0;i?) = iV 



(1 -I- cos^ 



R 



^tan-). 



where A'^ is a normalization constant given by 



N 



1 R^ 
R^' + - 



4 ' 



-1/2 



(6) 



(7) 



being ip^x/R) the mother wavelet associated to the conti- 
nous transform on the plane. R defines the size of the Mexi- 
can Hat. Then, the generalization to 52 of the MEXHAT on 
the plane is clearly given by eq. (6) with ■<p{x/R) given by 
eq. (4). Therefore, moving any point to the North Pole and 
taking polar coordinates, the wavelet coefficients associated 
to the function f{9, 0) are given at that point by 
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Figure 1. Analysed co-added COBE-DMR map (top) and the corresponding wavelet coefficient map obtained at 4.6 degrees (bottom). 
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Figure 2. Skewness values obtained from the analysed co-added COBE-DMR map at several scales are represented by dots in the top 
figure. The bands represent the 68%, 95% and 99% confidence levels. The corresponding kurtosis results are presented in the bottom 
figure. Values at zero correspond to real space. 



w{R) = / dt7/(6»,0)*s2(6';i?). (8) 



3 METHOD 

The methodology used in this work is designed to anal- 
yse the 4-year COBE-DMR data in the HEALPix pixelisa- 
tion. Maps at 53 and 90 GHz are added using inverse-noise- 
variance weights. The resolution of these maps is Nside = 64 
corresponding to a pixel size of « 55' and resulting in 49152 
pixels. Only those pixels surviving the extended Galactic cut 
of Banday et al. (1997) are considered (in order to excise 
those pixels with significant Galactic contamination near 
the plane of the Galaxy). Moreover, we subtract the best- 
fit monopole and dipole (as computed on the cut-sky) from 
the remaining pixelsQ We refer to the resulting map as the 
co-added COBE-DMR map. 

The spherical Mexican Hat wavelet of several sizes is 

t Note that this differs slightly from the treatment by Barreiro et 
al. It has been checked that in the formalism used in that paper, 
removing monopole and dipole does not affect the results 



convolved with the co-added COBE-DMR map. The size 
of the Mexican Hat (R) fixes the characteristic scale of 
the wavelet coefficients (obtained after the convolution). 
Wavelet coefficient maps are obtained for several Mexican 
Hat sizes and statistical tests presented below are applied 
to these maps (see Figure 1 for an example of a wavelet 
coefficients map). Confidence intervals are estimated from 
simulations. We simulate CMB temperature fiuctuations as- 
suming a CDM power spectrum (computed using the CMB- 
FAST software) with D,t = 0.05, D.c = 0.3, flv = 0.65, 
ill/ = 0, Ho = 65 km/s/Mpc and n = 1. The power spectrum 
is normalized to Qrms-ps = ISf^K. Simulated maps are 
smoothed with a Gaussian beam of FWHM = 7.032° (to 
approximate the beam response of the DMR instrument). 
Noise maps, generated with the same observation pattern 
and rms noise levels as the COBE-DMR data, are added to 
the simulated signal. The Galactic mask is also applied to 
the simulations, and the monopole and dipole components 
are estimated and subtracted outside the excluded region 
for each simulated map. Since the Mexican Hat wavelet ex- 
tends to the whole sphere, convolution with the sky at a 
fixed pixel will involve information from pixels over all an- 
gular separations. However, in this analysis the wavelet coef- 
ficients must be computed excluding those pixels within the 
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Galactic mask. In principle, this effect could bias the distri- 
bution we want to test. In practise, there is no appreciable 
change in results introduced by the partial sky coverage. 

We have applied several tests in wavelet space that can 
be divided into two categories. The first ones are computed 
globally from wavelet coefficient maps at scales ranging from 
55' to 1760'. We have also computed these tests on the real 
map and results are indicated by scale 0' in Figure 2 and 
Table 1. In addition we have also considered tests based 
on the derived distribution function at each pixel, at scales 
ranging from 108' to 1728'. In this way we take into account 
the different amplitude of the noise at each pixel. The rea- 
son for considering scales different from the ones used in the 
previous case, is that the maximum resolution of the anal- 
ysed COBE-DMR map and of the simulations is Nsizn = 32 
(corresponding to 12228 pixels on the sphere with size 1°.8). 
Information on temperature distributions at each pixel has 
to be kept at each scale. We needed to degrade the resolution 
to be able to handle all this information. 

In the first category of tests we calculated different 
statistics commonly used to detect non-Gaussianity. Skew- 
ness and kurtosis will test the third (ks) and fourth (fC4) 
order cumulants : 



S ■ 



I 3/2 
US,/ 1^2 



K ■ 



K4/K2 



(9) 



Correlations between two consecutive scales k,j are tested 
using the conventional 2-scale correlation given by: 



J2^^^AT/T{k,i)AT/T{j,t) 
Nr,V[k)V[j) 



where 



V{k) = ( 



Et^i(Ar/r(fc,i))^ 



,1/2 



(10) 



(11) 



and Np will correspond to the number of pixels considered at 
the largest scale (the ones outside the Galactic mask located 
at a distance from the border larger than twice the Mexican 
Hat size corresponding to that scale) . As well as a correlation 
analogous to the one used by Barreiro et al. 2000, given by: 



C2,2(fc,j) = 



Y,'^:^{AT/nk,i))\AT/T{j,i)f 



(12) 



V{k)V{j) 

We have also considered correlations between three consec- 
utive scales k,j,l given by: 



C^{k,j,l) 



j^N,^ Ar/r(fc, i)AT/T{j, i)AT/T{l, i) 
NpV{k)V{j)V{l) 



(13) 



In addition we compare the temperature fluctuation distri- 
butions from the data to those obtained from simulations, 
at several scales. The Kolmogorov-Smirnov test is used in 
this comparison. We have calculated the distance parameter 
for the COBE-DMR data as well as a distance parameter 
distribution from simulations. 

In the second category of tests we have calculated 
wavelet coefficient distributions at each pixel at each scale. 
We compute the probability of extrema taking into account 
the dispersion at each pixel. The second of these tests con- 
sists in the following. At a fixed scale, for the COBE-DMR 
data we select those pixels with values outside the 99% and 
95%. The number of those selected pixels is compared with 
the expected number from simulations. 



4 RESULTS AND CONCLUSIONS 

We have applied the statistical tests described in the pre- 
vious section to the COBE-DMR data. 500 Gaussian sim- 
ulations were performed to obtain the confidence intervals 
for each test as well as the distribution functions at each 
pixel. The results are summarised in figures 2 and 3 and 
tables 1, 2 and 3. For the tests in the first category, cor- 
responding to figures 2 and 3, all but one of the COBE- 
DMR values are inside the 95% confidence interval. In the 
case of C2,i between scales 55' and 110', the value for the 
COBE-DMR is inside the 99% confidence interval. Since the 
correlations are calculated between scales that do not differ 
from each other by a constant amount, one can observe an 
oscillating behaviour for C'2,1 and (72,2- Some peaks occur 
when the correlation is calculated between two scales dif- 
fering by a small amount. Also part of the first category 
tests, the Kolmogorov-Smirnov distance obtained from the 
COBE-DMR data is, at all scales, smaller than the maxi- 
mum found for the 500 Gaussian simulations (Table 1). In 
the second category of tests (Tables 2 and 3), only the prob- 
ability of finding a number of pixels outside the 99%, 95% 
confidence intervals, greater or equal the number of COBE- 
DMR pixels found outside those intervals, at the smallest 
scale (108'), is very low. However, noise might be dominat- 
ing at this scale. Therefore, no obvious deviation from Gaus- 
sianity can be noted from any of the tests presented in this 
work. 

The spherical Mexican Hat wavelet is introduced in this 
paper as a tool to study the information at different scales 
in the COBE-DMR data. The only other wavelet on the 
sphere previously used to analyse CMB data is the spheri- 
cal Haar wavelet (Tenorio et al. 1999, Barreiro et al. 2000). 
The main difference between these two wavelets is the geo- 
metrical property of isotropy that characterises the Mexican 
Hat. The Haar wavelet however biases certain directions. 
As a consequence, orientation of the data has to be taken 
into account in case this last wavelet is used in the analy- 
sis. The formalism of projecting the Mexican Hat (defined 
on the plane) on the surface of the sphere was presented 
in Section 2. Specific details related for example to the im- 
plementation on the HEALPix pixelitation will be given in 
Martinez-Gonzalez et al. (2001). 

In summary, several statistical tests have been per- 
formed in wavelet space. They are divided in two cate- 
gories. The first one includes tests performed globally on 
the wavelet coefficient maps. Tests in the second category 
take into account local properties at each scale. Tempera- 
ture fluctuation distributions are computed from Gaussian 
simulations at each pixel at each scale. The local dispersion 
is considered in the statistical tests applied. These methods 
have not been applied in any of the previous works. From 
all of them we can conclude that the analysed COBE-DMR 
data are compatible with Gaussianity. Skewness, kurtosis 
and C2,2 were also tested by Barreiro et al. 2000 with the 
spherical Haar wavelet. There is no indication from this work 
that one of the two wavelets is more optimal for detecting 
non-Gaussian features in the COBE-DMR data. However, 
the Mexican Hat wavelet is more localized in wavelet space 
than the Haar wavelet. Moreover, a comparison of the per- 
formance on the sphere of the two wavelets shows, that the 
spherical Mexican Hat wavelet is more optimal than the 
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Figure 3. Values corresponding to tlie 2-scale correlation C2,i are presented in the top figure. At each scale, the correlation has been 
calculated with the scale directly above (i.e. the value at 55' corresponds to the correlation between wavelet coefficients at 55' and 110'). 
The scales are the same ones appearing in Table 1, from 55' to 1100'. Bands correspond to confidence levels at 68%, 95% and 99%. The 
middle figure presents the results for the 2-scaIe correlation defined as in Barreiro et al. 2000 6*2,2. The figure in the bottom corresponds 
to the correlation between three consecutive scales. At each scale the value presented is the correlation between that scale and the other 
two direclty above that one (i.e. at 55' the correlation is calculated between wavelet coefficients at that scale and 110' and 220'). 



Table 1. Kolmogorov-Smirnov distance dKS 



scale(') 


dKS(COBE) 


Max dKS 





0.028 


0.035 


55 


0.036 


0.043 


110 


0.034 


0.041 


220 


0.033 


0.052 


275 


0.029 


0.068 


440 


0.036 


0.137 


550 


0.030 


0.173 


825 


0.088 


0.296 


880 


0.107 


0.351 


1100 


0.220 


0.560 
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Table 2. Probability of Extrema 



scale (') Vroh{T^Traa:.COBE)% Prob(T ^ T„,„COSE)% Vrob(\{T /up,^)\ ^ \{T /up,^)^axCOBE\)% 



108 84.2 80.6 96.8 

324 38.4 86.6 82.6 

432 15.4 90.6 41.6 

648 24.0 84.4 44.6 

864 67.6 62.4 78.2 



Table 3. Probability of Areas outside the 99% and the 95% 



scale (') COBE pix out 99% Prob(pix out 99% ^ COBE) % COBE pix out 95% Prob(pix out 95% ^ COBE) % 



108 87 0.0 359 2.0 

324 44 37.0 187 82.4 

432 34 48.4 161 77.6 

648 18 45.6 50 89.2 

864 100.0 22 82.0 
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